Noname manuscript No. 

(will be inserted by the editor) 



A new method of description of three-particle Coulombic 
systems 

V. B. Belyaev • A. A. Naumkin 



Received: date / Accepted: date 

Abstract We present a method for treatment of three charged particles. The pro- 
posed method has universal character and is applicable both for bound and continuum 
states. A finite rank approximation is used for Coulomb potential in three-body system 
Hamiltonian, that results in a system of one-dimensional coupled integral equations. 
Preliminary numerical results for three-body atomic and molecular systems like H~ , 
He, ppfj, and other are presented. 

Keywords three-body atomic systems ■ hyperspherical functions ■ finite rank 
operator ■ bound states 

PACS 03.65.Ge ■ 31.15.-p ■ 31.15.xj ■ 36.10.Ee 
1 Introduction 

The quantum three-body problem emerges in various fields of physics, and different 
methods of treating it are developed. However, there are no universal methods able to 
solve it in case of charged particles in the continuum. These problem is important in 
atomic and molecular physics and in nuclear astrophysics. 

The main purpose of the work is to develop a procedure applicable to treatment 
of three charged particles in continuum. One of the ways to construct it is to make an 
approximation on operator level, i.e. in the Hamiltonian of three-body system under 
consideration. After that any boundary conditions can be used. The idea of our work 
was inspired by earlier paper of A. Weinstein [1], who introduced the so-called method 
of intermediate Hamiltonians. This method was applied to calculate lower bounds for 
eigenvalues of some differential operators. N. W. Bazley and D. W. Fox applied it to 
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He atom and other physical systems [2][3]- They constructed sequence of intermediate 
Hamiltonians using finite rank operators. These operators are defined in the whole 
space of full Hamiltonian. In opposite to that we will apply finite rank approximation 
in a subspace of the three-body system Hamiltonian, namely in the angular space of 
hyperspherical variables. This results in a system of coupled one-dimensional integral 
equations. 

In the following sections we review a method of intermediate Hamiltonians, derive 
system of integral equations and report the results of calculations. 



2 Hyperspherical coordinates 

There are different ways to formulate the three-body problem. We use hyperspheri- 
cal coordinates, and in this section we give a brief review of them. Complete theory, 
derivations etc can be found, e.g., in [B]. 

We start with the three-body system Hamiltonian: 

3 i 

i— 1 i<j 

where i enumerates different particles and corresponding sets od Jackobi coordinates. is 
a numbr,; is a position vector of the i-th particle. The scaled Jackobi coordinates are 
introduced as follows: 
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mi + 77i2 + JTI3 
and the Hamiltonian fTJ takes the form 



mfrj + m k r k 
raj + m k 



(3) 



where V = V13 + V23 + V31 — a sum of pair potentials. Taking x and y in spherical 
coordinates (x,y) — >• (x,9i, fi,y, O2, 922)1 one obtains: 
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Now let us introduce hyperspherical variables: 

x = pcosa, y — psina 



'■JL 

dx 
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Here p is hyperradius, a — hyperangle. Hamiltonian expressed in terms of this variables 
has the form: 
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Angular part of a kinetic energy operator is the hypermomentum operator: 

K = + 4 cot 2a^- + —L-A ni + —L-A n2 , (7) 
dor aa cos^ a sm a 



and its eigenfunctions are hyperspherical harmonics 

'■' h+ *\coB2a)Y hmi (fh)Y hma 

(8) 



where 

1/2 

(9) 



„M 2 



2n\(K + 2)(n + Zj + Z 2 + 1)! 



r(n + Zi + 3/2)r(n + h + 3/2) 



Let us consider system of three particles with masses mi,m2, 772,3 and charges Qi, Q2, 93- 
The Coulomb potential has the form: 

V(x,y) = ^ + ^ + ^, (10) 

XI X 2 X 3 



where bi = ./ \ _~ OjOfc. In hyperspherical coordinates: 

npJfi) = If^_+^+^) (11) 

p ycOSai COSQ2 COSQ3 / 

Here cti — hyperangles corresponding to different sets of Jacobi coordinates. 



3 Finite rank operators 

Finite rank operators are widely used in different problems of mathematical physics. 
They allow one to reduce complexity of a problem and proceed to its solution. E.g., 
in [5] finite-rank operator was used to describe nuclear part of full Hamiltonian in a 
problem of low energy it — 3 He scattering. 

N. W. Bazley and D. W. Fox used finite rank operators to calculate lower bounds 
of eigenvalues of Schrodinger eqution 2,3 . Let us shortly review the method of inter- 
mediate Hamiltonians they used. 

We suppose that full Hamiltonian H can be presented as a sum of H°, that has 
known eigenvalues and eigenfunctions, and a positively definite H' . The exactly solvable 
Hamiltonian H is assumed to have ordered discrete energy levels El < < ... below 
its continuum spectrum. The corresponding eigenfunctions are i/i®, and we have 

H°iP? = EUl (12) 

Since H = H° + H' , where H' is positively definite, H° < H and < E\. Thus, the 
full Hamiltonian H and H are linked by a sequence of intermediate Hamiltonians: 

H° < H k < H k+1 < H. (13) 

To construct the Hamiltonians H k , we introduce a system of k linearly independent 
functions Pi,P2, The set of functions pi,P2, ■■■ is defined in the whole space of 
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definition of the Hamiltonian H. Projection of some wavefunction ip on these functions 
is given by 



P k ip = ^a kPk (14) 



The projection P k increases with k: 



< <¥>|P>) < (p\P k+ \) < Mp) (15) 
< (p\H'P k p) < (p\H'P k+1 p) < {tp\H'tp) (16) 

From Eq. (JTSJ) we can see that H'P k < H'P k+1 < H', and we now define intermediate 
Hamiltonian as 

H k = H° + H'P k (17) 

It is important to to emphasize that the finite rank operator H k acts on the functions 
Pi,P2, •■■ in the same way as full Hamiltonian H: 

H k \i)=H\i), i = l,...,k (18) 

This is the main property of a some finite rank operators which we will use below. 

Following this idea, we construct such an operator in the angular space of defini- 
tion of operator (0. The Coulomb potential in hyperspherical variables has the form 
V(p,D) = where f(f2) is the angular part of potential. We use a finite rank 

approximation in it. Namely, the function f((2) is replaced by a finite rank operator: 

N 

f(n)^f N =Y,f\<Pi)d ij (<P j \f (19) 
*>J 

Here <pj are some auxiliary functions defined in angular space, dy = ((p^f^j)^ 1 - 
inverse matrix element. 



4 Formalism 

Here we derive a system of coupled one-dimensional integral equations using the finite 
rank approximation. Let us start with the Schrodinger equation: 

(H + V)\&)=E\P) (20) 

Here H$ is kinetic energy, V — interaction potential. For bound states this equation 
can be written in integral form using free Green function: 

\W) = (E-H y 1 V\&) = -G E V\V) (21) 

Let us rewrite it in coordinate representation: 

<P(R) = -J dR'G B (R,R')V(R')* r (R'), (22) 

where R = (x,y) = (p,fl), and use the Coulombic potential: V(R) = i/(X2). We 
obtain integral equation for the wavefunction in hyperspherical coordinates: 

9{j>,0) = - I / ''dp 'dSl 'G E (p,p '■Q,Q')^f{Q')^(p ',fl') (23) 
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Using finite rank operator (|19[) instead angular part of potential f(O), we obtain rep- 
resentation for the wavefunction <P": 



N 



V{p,n) = -J2J P dp'dtfGjsfap'; n,ff)f dijCjQ/), (24) 

where Cj(p') = J dtt" <p j (tt")f(n")t'(p, n") — new unknown functions. 

In order to obtain a system of integral equations for functions Ci(p), we use integral 
operator: J dfl ip k (ft)f(n) .... As a result we obtain: 



C k(p) = -Y,j d P P I dndn'^ k {Q)f{Q)G E { P ,p' ] Q,n')f{Q') i p i {Q')d ij C j {p') 

(25) 

or 

G k (p) = -J2 I dp' M kl (p, p')dij Cj(p'), (26) 
M kl (p,p') = p' 4 I ' dQdQ' ' <p k (n)f(I2)G E (p,p';n,lf (27) 



The Green function G E (R-, R-') has the simplest form in the momentum represen- 
tation. Using the plane wave expansion: 



e iqx+w = 7-7T2 E i K JK + 2^ P )y^i M {n p )y% l i M {n K 



(2tt) 3 (k P ) 



dpdq . 2m /h 
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G_b(x, y) = / / -— g exp(ipx + iqy 



yy (27r)° p z + q z + k z 

we obtain partial harmonics of free Green function: 

g% ( P , P ') = J J yk l l M (o) G E (K, R') y l ^ l l M (Q K )dQ dn' 
= 1 ^(p) -WM-W/V)^ 



2m£ 



/°'\ J ^K+2(«:op)-R'k+2(^op'). < p < p 



(2tt)3 V p ) {K k+2 (kop)I k +2(kop'), < p' < p 

Here J n (x), In{x) and K n (x) are Bessel function and modified Bessel functions of first 
and second kin< 
equations (l25"Tl : 



and second kind, respectively, k 2 = p 2 + q 2 . Now we can calculate the kernels of integral 



M ki (p, P ')= G%(p,p')^ k \f\y l ^i M ){y^i M \f\n) (28) 

KLMhh 

We derived a system of coupled one-dimensional integral equations (|26|l . Now one need 
to calculate the kernel and solve this system numerically. At this stage of treating the 
Coulomb three-body problem the finite rank approximation makes it sufficiently easier. 
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5 Calculation and results 

We constructed the finite rank operator (|19[) using hyperspherical functions. They have 
been chosen for convenience, but one can use some other set of linearly independent 
functions defined in angular space. 

It is important to mention that the representation (|25p for solution of Schrodinger 
equation is not a well known hyperspherical expansion. One can see it from the defini- 
tion of Ci (p) . 

We performed calculations using finite rank operators constructed on 1, 3 and 

6 auxiliary functions. In calculation the kernel (|28p we should summate an infinite 
number of terms, but we stopped at values of the hypermomentum K equal to 6, 10 
and 14. In order to solve integral equations, the variables p and p' were discretized 
with 100 mesh points. 

We calculated binding energies of the ground state of such systems: He, H~ , 
ppp and ddp. Results of these calculations are presented in Table [1] 



Table 1 Calculated and precise binding energies, eV 





E pr , eV 


FRA 


Kmax — 6 


^max — 10 


i^max — 14 






N=l 


25.0 


20.0 


18.2 




14.34 


N=3 


20.0 


18.5 


17.1 






N=6 


18.0 


16.2 


15.6 






N=l 


102 


98 


95 


He 


79.0 


N=3 


99 


91 


89 






N=6 


95 


87 


85 






N=l 


7.3 


10.0 


11.0 


vi 


16.25 


N=3 


8.5 


12.0 


13.7 






N=6 


10.1 


13.5 


15.1 






N=l 


1050 


1700 


1850 


ppp 


2782 


N=3 


1360 


2044 


2101 






N=6 


1690 


2290 


2332 






N=l 


1200 


1820 


1990 


ddp 


2988 


N=3 


1540 


2072 


2480 






N=6 


1845 


2195 


2654 



The precise energies are taken from [7]. 



6 Conclusion 

Binding energies of different three-body Coulombic systems were calculated within 
a finite rank approximation method. The finite rank approximation is made in an 
angular part of potential in three-body Hamiltonian. This method was tested on some 
of these systems earlier in [4]. The results obtained shows it can be useful for solving 
the Coulombic three-body problem. 

One can expect that this method will also be applicable to three charged particles 
in continuum, since the approximation (|19p is made at the operator level, but not in 
the wavefunction. 
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